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ABSTRACT 


This  report  deals  with  the  geometrical  theory  of  diffraction  (GTD) 
solution  to  the  problem  of  plane  wave  high  frequency  (h.f.)  scattering  by 
a  smooth,  perfectly-conducting  cylindrical  surface.  According  to  the  GTD, 
the  surface  diffracted  rays  accurately  describe  the  field  behavior  deep  in 
the  shadow  region  behind  the  obstacle;  whereas,  the  incident  and  reflected 
rays  adequately  describe  the  field  in  the  lit  region  sufficiently  far  from 
the  optical  shadow  boundary.  These  ray  solutions  fail  within  the  trans¬ 
ition  region  adjacent  to  the  shadow  boundary;  hence,  an  approximate  h.f. 
solution  is  developed  for  describing  the  field  within  the  transition  region. 
The  usefulness  and  the  accuracy  of  this  transition  region  solution  is 
illustrated  by  analyzing  the  h.f.  plane  wave  scattering  by  a  2-D  cylindrically 
capped  thick  half  plane,  and  a  2-D  circular  cylinder  (for  which  an  exact 
solution  is  available  for  comparison),  respectively.  The  present  transition 
region  solution  is  in  a  form  which  facilitates  its  generalization  to  treat 
2-D,  and  3-D  convex  surfaces  of  non-constant  curvature.  A  simple  example 
of  this  generalization  is  illustrated  by  analyzing  the  h.f.  plane  wave 
scattering  by  a  3-D,  hemi spherically  capped  semi -infinite  circular  cylinder. 
The  work  described  in  this  report  is  important  to  the  analysis  of  the  problem 
of  h.f.  scattering  by  a  finite  elliptic  cylinder  excited  by  a  finite  source; 
this  problem  is  directly  relevant  to  a  study  of  the  effect  of  an  aircraft 
fuselage  on  the  performance  of  wing  mounted  antennas. 
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I. 


INTRODUCTION 


On-aircraft  antenna  pattern  performance  is  very  much  dependent  on 
the  scattered  field  from  the  fuselage.  This  point  was  emphasized  in  the 
analysis  of  fuselage  mounted  antennas[l].  For  certain  sectors  of  the 
volumetric  pattern,  the  fuselage  may  even  play  a  more  dominant  role  for 
antennas  mounted  off  the  fuselage.  This  is  especially  true  in  the 
azimuth  plane  for  wing  mounted  antennas.  For  many  years  scale  model 
measurements  have  been  used  to  study  these  effects;  however,  such 
measurements  are  both  time  consuming  and  expensive.  On  the  other  hand, 
numerical  solutions  have  been  successfully  applied  for  the  past  few 
years  to  study  designs  and  locations  for  fuselage-mounted  antennas. 
Fundamentally,  the  analyses  developed  for  fuselage-mounted  antennas  can 
be  extended  to  handle  off- fusel  age  antennas. 

Based  on  the  numerical  techniques  previously  developed[l  ,2,3],  it 
is  proposed  that  similar  solutions  be  investigated  for  wing-mounted 
antennas.  Specifically,  the  near-zone  scattered  field  solutions  for  a 
finite  flat  plate  and  finite  circular  cylinder  have  already  been 
developed[l ] .  However,  these  models  are  not  as  versatile  as  desired. 
Further  the  finite  cylinder  is  limited  to  electrically  small  radii  due 
to  use  of  a  modal  solution.  Consequently,  it  is  proposed  that  a  high 
frequency  solution  for  sources  in  the  near  zone  of  a  finite  elliptic 
cylinder  be  employed.  This  model  allows  for  variations  in  the  cross- 
section  of  the  theoretical  model  and  is  not  limited  to  lower  frequencies, 
i.e.,  to  electrically  small  bodies. 

The  theoretical  solutions  necessary  to  analyze  this  problem  are 
developed  in  this  report  in  terms  of  the  solution  to  the  scattering  by  a 
two-dimensional  circular  cylinder.  The  geometric  optics  solution 
satisfactorily  handles  the  lit  zone;  whereas,  the  geometrical  theory  of 
diffraction  can  be  used  in  the  deep  shadow.  However,  both  solutions  fail 
on  and  near  the  shadow  boundary;  hence,  an  approximate  asymptotic  high 


frequency  result  is  developed  for  analyzing  the  field  within  the  transition 
region  adjacent  to  the  shadow  boundary.  Further,  the  solutions  developed 
herein  can  be  extended  to  handle  arbitrary  curved  surfaces  including  the 
three-dimensional  case.  These  solutions  will  provide  the  necessary 
theoretical  background  to  complete  the  study  of  a  finite  elliptic  cylinder 
which  will  be  treated  in  a  future  report. 

The  canonical  problem  that  is  treated  involves  the  calculation  of 
the  total  high  frequency  field  which  is  observed  at  a  finite  distance 
from  a  convex,  cylindrically  shaped  diffracting  obstacle  when  it  is 
illuminated  by  an  incident  plane  wave.  This  problem  is  directly 
equivalent  to  the  relevant  problem  of  calculating  the  radiation  pattern 
of  a  source  located  at  a  finite  distance  from  the  same  smooth,  convex 
diffracting  obstacle  via  the  reciprocity  theorem  for  electromagnetic 
fields. 

II.  HIGH  FREQUENCY  SCATTERING  BY  SOME  SPECIAL 

2-D,  SMOOTH,  CYLINDRICAL  SURFACES 

This  section  deals  with  the  high  frequency  scattering  of  electro¬ 
magnetic  waves  by  a  smooth,  two-dimensional  perfectly-conducting  cylindrical 
surface.  In  this  study,  special  emphasis  is  given  to  the  region  near  the 
shadow  boundary.  A  solution  is  desired  which  can  be  generalized  (based  on 
high  frequency  approximations)  to  treat  the  scattering  by  three-dimensional 
surfaces  without  having  to  solve  additional  canonical  problems.  The  general 
problem  of  the  diffraction  of  waves  by  smooth,  convex  surfaces  for  grazing 
angles,  (i.e.,  in  the  transition  region  adjacent  to  a  shadow  boundary)  is 
difficult  to  solve  and  it  has  been  the  subject  of  investigation  by  various 
authors[4]  who  considered  some  general  cases;  however,  owing  to  the  mathe¬ 
matical  complexity  of  their  solutions,  their  results  do  not  appear  to  be  in 
a  form  tractable  for  application  to  engineering  problems  of  practical 
interest.  For  the  special  case  of  high  frequency  plane  wave  diffraction  by 
a  perfectly-conducting  circular  cylinder,  Goriainov[5]  has  obtained  results 


for  grazing  angles  which  are  in  a  form  suitable  for  numerical  calculations. 
Wait  and  Conda[ 6]  extended  Goriainov's  high  frequency  results  to  treat  the 
plane  wave  diffraction  by  a  circular  cylinder  with  an  impedance  type 
boundary  condition;  they  also  included  the  case  of  spherical  wave  incidence. 
Most  of  the  previous  works  on  curved  surface  diffraction  for  grazing  angles 
including  the  work  in  References  [5]  and  [6]  make  use  of  Fock  integrals 
introduced  by  Fock  in  his  fundamental  work  on  diffraction  by  smooth  surfaces 
(many  of  Fock's  papers  are  listed  in  Reference  [4]).  Wait  and  Conda[6] 
indicate  that  their  high  frequency  (h.f.)  solution  for  the  plane  wave 
diffraction  (for  grazing  angles)  by  a  circular  cylinder  can  be  interpreted 
as  being  the  superposition  of  the  result  based  on  Kirchhoff  diffraction 
theory,  and  a  second  term  which  is  viewed  as  a  correction  to  the  Kirchhoff 
theory.  This  observation  by  Wait  and  Conda[6]  constitutes  the  physical 
basis  for  the  development  pursued  in  this  report  to  arrive  at  a  result 
which  could  be  simply  extended  to  treat  the  h.f.  scattering  (in  the 
transition  zone)  by  arbitrary  2-D  and  3-D  smooth,  convex,  conducting  surfaces. 
The  decomposition  of  the  field  scattered  by  the  2-D  cylinder  into  a  sum  of 
a  Kirchhoff  diffraction  term  and  a  surface  diffracted  term  (viewed  as  a 
correction  to  Kirchhoff  theory)  as  observed  by  Wait  and  Conda[6]  was  also 
indicated  by  Fock  (see  pp.  194-196  of  Reference  [6]).  In  his  work  on  the 
3-D  scattering  problem  of  the  Fresnel  diffraction  by  a  sphere  for  the  special 
case  of  source  and  observer  raised  to  small  heights  above  the  sphere,  Fock 
expressed  the  scattered  field  in  terms  of  a  superposition  of  two  integrals; 
one  of  these  integrals  happened  to  be  independent  of  the  electrical  properties 
of  the  surface,  and  the  other  explicitly  contained  the  surface  properties 
and  its  associated  electrical  properties.  Fock  then  showed  that  the  integral 
for  the  scattered  field  which  is  independent  of  the  electrical  properties 
reduced  to  a  Fresnel  edge  diffraction  term;  this  term  corresponds  to 
Kirchhoff  diffraction.  On  the  other  hand,  the  integral  which  depended  on 
the  electrical  properties  of  the  surface  was  shown  by  Fock  to  yield  the 
surface  diffracted  term  (which  is  viewed  as  a  correction  to  Kirchhoff 
diffraction).  Consequently,  it  is  believed  that  the  basis  of  our  analytical 
development  given  in  this  report  is  justified  in  view  of  the  above  remarks. 


Since  the  circular  cylinder  has  two  penumbral  regions  on  its  surface 
as  shown  in  Fig.  1,  one  needs  to  isolate  the  contribution  to  the  field 


Fig.  1.  Geometry  for  the  plane  wave  scattering 
by  a  circular  cylinder. 

from  each  of  these  penumbral  regions  in  order  to  allow  an  extension  of 
the  circular  cylinder  analysis  to  treat  the  h.f.  waves  diffracted  from 
a  single  penumbral  region,  such  as  that  which  occurs  on  a  smooth  convex 
surface  of  the  nose  of  an  aircraft  when  it  is  illuminated  by  a  wing 
mounted  antenna.  One  notes  that  the  angular  extent  of  the  penumbral 
regions  and  the  transition  regions  (which  lie  adjacent  to  the  shadow 
boundaries)  is  approximately  2(ka/2)~^  in  radians  (ka  =  circumference 
of  the  cylinder  in  wavelengths).  In  Section  I I -A  below,  it  is  shown  that 
the  h.f.  field  contributions  arising  from  each  of  the  two  penumbral  regions 
on  a  circular  cylinder  can  be  identified,  and  separated  in  a  manner  which 
enables  one  to  calculate  the  h.f.  field  within  a  transition  region  which  is 


4 


associated  with  a  single  penumbral  region  on  the  cylinder.  Such  a  decompo¬ 
sition  of  the  circular  cylinder  solution  facilitates  the  analysis  of  the 
problem  of  h.f.  plane  wave  diffraction  by  the  convex  portion  of  the  2-D, 
perfectly-conducting,  rounded-end,  thick-screen  configuration  of  Fig.  2. 


TRANSITION 


Fig.  2.  Geometry  associated  with  the  plane  wave  scattering 
by  a  cylindrically  tipped  thick  half-plane. 

One  notes  that  the  configuration  of  Fig.  2  consists  of  a  half  cylinder  of 
radius  (a)  which  is  smoothly  joined  to  a  thick  half  plane  whose  thickness 
is  2a.  The  curve  defining  the  boundary  of  the  configuration  in  Fig.  2  is 
not  closed,  and  therefore  there  exists  only  one  penumbral  region  on  the 
convex  portion  of  this  scattering  geometry  as  illustrated  in  Fig.  2.  The 
results  for  the  h.f.  scattering  by  the  canonical  configuration  of  Fig.  2 
are  given  in  Section  1 1 -B ;  it  is  conjectured  that  these  results  can  be 
extended  in  a  systematic  fashion  to  treat  the  h.f.  scattering  by  arbitrary 
convex  conducting  surfaces  if  one  employs  the  hypothesis  valid  at  high 
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frequencies,  namely  that  h.f.  diffraction  and  reflection  are  local 
phenomena  (i.e.,  h.f.  diffraction  and  reflection  mechanisms  are  localized 
to  the  immediate  vicinity  of  certain  points  on  the  surface  of  the  scatterer, 
which  correspond  to  points  of  diffraction  and  reflection,  respectively).  An 
extension  of  the  results  in  Section  II-B  for  the  2-D  canonical  configuration 
of  Fig.  2  to  analyze  the  problem  of  h.f.  scattering  by  a  3-D,  perfectly- 
conducting  hemi spherically  capped,  semi -infinite  circular  cylinder  geometry 
is  indicated  in  Section  III.  The  more  general  extension  of  the  results  of 
Section  II-B  to  handle  arbitrary  curved  surfaces  will  be  the  topic  of  a 
future  investigation. 

A.  Transition  Region  Solution  for  the  H.F,  Plane  Wave 
Scattering  by  a  Circular  Cylinder 

Let  the  total  field,  u  in  a  region  external  to  the  2-D  perfectly- 

conducting  circular  cylinder  of  Fig.  1  be  the  superposition  of  the  incident 

field  (u  )  and  the  scattered  field  (u  ).  In  the  case  of  TMz  excitation, 

the  electric  field  associated  with  the  incident  plane  wave  (u1 )  is  entirely 

parallel  to  the  axis  of  the  cylinder  (i.e.,  along  z);  furthermore,  the 

total  electric  field  denoted  by  u  is  also  entirely  z -directed.  Similarly, 

for  the  TEz  excitation  of  the  cylinder,  u1  and  u  represent  the  incident 

and  total  magnetic  field  intensities,  respectively,  which  are  entirely 

z-directed.  Let  u„  denote  the  value  of  u  for  the  TM,  case,  and  let  u. 

s  L .  h 

denote  the  value  of  u  for  the  TEz  case.  Then,  us  =  u^+s  satisfies  the  2-D 

h  h 

homogeneous,  reduced  wave  equation  in  the  region  exterior  to  the  cylinder, 
i.e., 

(1)  (vj  +  k2)  us  =  0, 

h 

subject  to  the  boundary  conditions; 

(2a)  us  =0  on  the  cylinder, 

or 
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0  on  the  cylinder. 


In  addition. 


(3)  tig  must  satisfy  the  Sommerfeld  radiation  condition  as  p-*-*  for  an 

h 

eiwt  time  convention  which  is  assumed  and  suppressed. 

The  incident  plane  wave  field  is  given  by: 

(4)  uj  =  eikx. 

h 

The  quantity  k  present  in  the  above  equations  refers  to  the  free  space 
wave  number.  Since  Eq.  (4)  satisfies  the  reduced  wave  equation,  it 

follows  from  Eq.  (1)  that  u|  must  also  satisfy  the  reduced  wave  equation. 

h  , 

An  application  of  the  2-D  Greens  theorem  to  the  region  surrounding  the 
cylinder  of  Fig.  1  yields 

...  <..  (2*r.  sg(pip')  ..  .  wD.il 


u(P)  •  u'(P)  +  J  u(P’)  JJ?- 


-  G0<plp,)  ^ 


where  G  (P|P')  is  the  2-D  free-space  Green's  function: 

(6)  Go(P!P')  =  -  j-H^2)(k]p  -  p'|). 

The  parameter  (p")  denotes  the  field  point  position  vector  of  P  =  P(p,<|>); 
and  p1  denotes  the  cylindrical  surface  position  vector  of  P‘  =  P 1  (p  ’,<}>' ) 
with  p'  =  a  on  the  cylinder.  Clearly,  the  integral  term  on  the  RHS  of 
Eq.  (5)  is  the  scattered  field  us.  Then,  from  Eqs.  (2),  (3)  and  (5) 


s  f2ir 

Ug(P)  =  -  Gp(P|P') 

°  J  0 


3US(P') 


ad4> ' , 


c  f2ir  3G„ 

uh<p)  ■  J0  k(p,) 


BG„(P|P') 


ad$  1 , 


with  the  understanding  that 


(8a)  us  =  us  +  us’  ^or  the  ™2  case* 


(8b)  u^  =  u^  +  u^,  for  the  TE^  case. 

It  is  well  known[4,5]  that  a  valid  asymototic  h.f.  approximation 
for  the  surface  field  which  is  induced  within  the  penumbral  regions  on 
the  surface  of  a  circular  cylinder,  when  it  is  illuminated  by  a  plane  wave, 
is  given  by  the  following  set  of  equations: 


3u  ~]  *  H.  .  -ika^T  .  .  -ikaiK 

— p  =a  *=u 


ika»p,  .  -ikaip0 
+  o(eJ  e 


WP'=a*Mri>e  *gUz) 


The  functions  which  appear  above  are  the  universal  Fock  functions  defined 
by 


(10a) 


^  -  If  P"i5t  -  1  e'iST 

» (10b)  g(t)  =  —  ]  iTTrT dx  »  9(0=  —  wttt)  dT 

A  Jr,  w2't  /T  Jr,  w2't 


in  which  the  fock-type  Airy  function  w2(t),  and  its  derivative  w2(t)  are 
defined  by: 


Ola), (11b) 


w2<t>  =  7:  1  e  T  dz  I  w2(t)  =  7:  |  z 
c  vT  Jr0  L  ^  ;r„ 


The  Fock  parameters  5-j  and  £2  are  given  by: 


(12a), (12b)  L  =  (|^)1/3  *.  ;  L  =  (l£)1/3  L  , 


8 


with  ij>i  =  j 4> 1  1  -  tt/2  +  2-ni,  and  ^  =  3tt/2  -  |<t>'  |  +  2itA.  The  contours  of 
integration  and  for  the  integrals  in  Eqs .  (10a, 10b)  and  (11a, lib) 
are  shown  in  Fig.  3. 


Fig.  3.  Contours  of  integration  for  the  Fock  and  Airy  functions. 

The  two  terms  on  the  RHS  of  Eq.  (9a)  and  Eq.  (9b)  are  commonly  viewed 
as  two  oppositely  traveling  creeping  waves  (associated  with  the  surface 
current)  on  the  cylinder  which  have  encircled  the  cylinder  t  times.  If 
(ka)  '  ir  »  1 ,  the  effect  of  the  multiply  encircling  terms  is  negligible 
and  the  inclusion  of  only  the  £=0  term  is  sufficient.  In  this  report,  one 
is  concerned  with  eventually  applying  the  appropriate  circular  cylinder 
results  to  the  semi-infinite  geometrical  configuration  of  Fig.  2;  hence, 
it  is  meaningful  to  retain  only  the  t=0  term  in  the  present  analysis. 


Employing  the  far-zone  approximation  to  6  (P  | P * )  leads  to  the 

3G  (PlP1)  0 

following  results  for  GQ ( P | P 1 )  and  - : 


03) 


[G„(P|P')1  *  -  i 


2  -i  ( kp -n/4 )  _i kacos  (<j> 1 ) 
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(14) 


3Go(pip,>1  „  k  nr 

3p‘  ,  4  lirkp 

p  -a 


cos  ($-<!> 1 ) 


e~i  (kp-ir/4) 


gikacos  ’ ) 


The  far-zone  approximation  implies  that  k|p  -  p'|  »  1  and  p  »  (p'=a)  in 
Eqs.  (7a, 7b). 


It  is  now  possible  to  obtain  asymptotic  approximations  to  the 
integrals  in  Eq.  (7a)  and  Eq.  (7b)  for  the  scattered  fields  subject  to 
the  restrictions  indicated  in  Eqs.  (13)  and  (14).  Employing  Eqs.  (9a) 
and  (13)  in  Eq.  (7a),  and  expanding  the  exponential  part  of  the  integrand 
about  a  stationary  point  allows  one  to  write  the  integral  in  terms  of 


functions  which  are  tabulated;  thus,  one  is  led  to  the  following  result  for 


A  similar  asymptotic  solution  for  the  integral  in  Eq.  (7b)  is  obtained  by 
incorporating  Eq.  (9b)  and  Eq.  (14)  in  it's  integrand;  thus. 


The  results  indicated  in  Eq$.  (15)  and  (16)  were  obtained  earlier  by 
Goriainov[5],  and  Wait  and  Conda[6].  The  universal  functions  p*U)  and 
q*U)  are  the  complex  conjugates  of  the  p(0  and  q(0  functions  defined 
by  Logan[4,5,7]  (p*(0  =  GU.q*-)  eiff/4  and  q*U)  =  GU,q=0)  el7r/4, 
where  G(4,q)  has  been  introduced  by  Wait[6]).  The  functions  e'^^4  p*(0. 


and  e~  '  q*($)  are  plotted  in  Figs.  4a  and  4b,  respectively 

arguments  C]  of  the  p*(5] )  and  q*(5-|)  functions  are  given  by: 


Fig.  4b.  Plot  of  e  q*(0  vs  t  based  on 
Logan's  tabulated  data[7]  for  q(£). 

(17a, 17b)  ^  =  (j£)1/3  ;  C2  =  (j£)1/3  i>2 

with  ip i  =  <j>  -  it  ,  and  'i>2  =  *  - 

The  details  of  the  asymptotic  approximation  to  the  integrals  in  Eqs.  (7a) 
and  (7b)  do  not  appear  to  be  readily  accessible  in  the  literature,  and  are 
therefore  given  in  Appendix  I  for  the  sake  of  completeness. 

It  has  been  pointed  out  by  Wait  and  Conda[6]  that  the  first  term  on 
the  RHS  of  Eqs.  (15)  and  (16)  corresponds  to  Kirchhoff  diffraction; 
whereas,  the  remaining  terms  may  be  viewed  as  a  correction  to  the  Kirchhoff 


theory.  The  Kirchhoff  result  is  independent  of  the  electrical  properties 

of  the  surface;  it  is  obtained  by  assuming  that  the  field  scattered 

by  the  cylinder  of  Fig.  1  is  equal  to  that  scattered  by  a  strip  of  the 

same  geometrical  cross-section  (i.e.,  strip  width  =  diameter  of  the 

cylinder)  which  is  located  at  x=0  as  in  Fig.  5a.  It  is  possible  to  provide 

a  ray  interpretation  for  the  scattered  fields  given  in  Eqs.  (15)  and  (16); 

the  Kirchhoff  term  is  composed  of  the  field  associated  with  the  Kirchhoff 

diffracted  rays  which  emanate  from  the  edges  and  (of  the  strip)  as 

illustrated  in  Fig.  5a;  whereas,  the  terms  involving  p*(£-|)  and  q*U-|) 

are  the  fields  associated  with  the  ray  which  is  launched  by  the  incident 

field  at  P-j  and  shed  along  the  forward  tangent  at  P^,  after  creeping  a 

distance  P^P3  on  the  cylinder  (P-j =  a  ^  | )  as  illustrated  in  Fig.  5b.  A 

pseudo-ray  interpretation^]  is  possible  for  the  terms  p*U2)  and  q*(^?) 

as  illustrated  in  Fig.  5b,  wherein  the  ray  associated  with  these  terms  is 

launched  at  ancl  is  shed  along  the  backward  tangent  at  P^  after  creeping 

a  distance  P^P^  on  the  cylinder  (P2P4  =  ak2l  =  a|^|);  this  ray  system 

unlike  the  ones  described  above,  does  not  satisfy  the  generalized  Fermat's 

principle  which  is  employed  in  Keller's  geometrical  theory  of  dif fraction[9] , 

but  it  is  nevertheless  useful  in  visualizing  the  geometrical  coordinates 

-ikaip„  -ikaiK 

which  are  associated  with  tne  terms  p*U9)  e  c  and  q*U9)  e  .  In 

k  ^  ^ 

Fig.  5a,  u,(P)  denotes  the  Kirchhoff  edge  diffracted  field  emanating  from 

k 

P-j ;  similarly,  u~(P)  denotes  the  field  arising  from  the  edge  P 2  so  that 

(18)  ufj(P)  +  Ug(P)  %  -  He  eW4  ~  .  ^r  *  $  * 

corresponds  to  the  Kirchhoff  result  which  is  contained  in  Eq.  (15)  and  Eq. 
(16),  respectively, 

One  is  now  able  to  identify  the  scattered  fields  which  arise  from 
the  two  isolated  penumbral  regions  shown  in  Fig.  1.  Specifically,  the 
scattered  field  at  P  (in  the  close  vicinity  of  ♦=*)  arising  from  the 
penumbral  region  associated  with  P^  is: 
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(19a)  u| (P ,P1 )  ^  u^(P)  -  e'11T/4  |X  (|^-)1/3  p*Uv 


-i k -j  e-ikP 


*"D 


for  the  TMz  case  and  <p  Z  tt 


MQkx  s,D  D  .  k,Dx  -i»/4  FT  /kaJ/3  »  'lkd,h  e'lkp 

(19b)  uh(P,P1)  ~  u-,  (P )  -  e  (jH  9  (€-,)  e 


/p 

for  the  TEZ  case  and  <p  %  v  . 


Similarly,  the  scattered  field  at  P  emanating  from  the  penumbral  region 
associated  with  is: 

(20a)  u|(P,Pa)  -  u|(P>  -  e-"'4  JT p*U2)  e',ka*2  o'"1' 


Si 


for  the  TM  case  and  <f>  ^  w 


(20b)  ^(P,P2)  -v  uk(P)  -  e'in/4{f  (|^)1/3  q*(52)  e  e 


/p 


for  the  TE_  case  and  4>  %  it  . 


Re-writing  the  RHS  of  Eq.  (18)  for  <j>-*Tr,  one  obtains  the  following: 

-ika  ( 4> -tt  ) 


(21) 


sinka(»-Tr)  ^(7~  iir/4  e_ikp 

(4>-ir )  virk 


(4>~tt  ) 


;ika(<t>-Tr) 
[i -it) 


\i_L  e'ilT/4  e 

Y2*k  e 


-ikp 


/p 


One  is  now  able  to  associate  the  first  term  on  the  RHS  of  Eq.  (21)  with 

.  J,, 

u.;  likewise,  the  second  term  on  the  RHS  of  Eq.  (21)  corresponds  to  iu. 

Ti  ka  (<f> -tt  ) 

However,  the  individual  terms  ±  - — —  on  °f  Eq.  (21)  are 

unbounded  at  <j>=ir  (even  though  the  LHS  of  Eq.  (21)  remains  finite); 
consequently,  they  must  be  replaced  by  the  bounded  Kirchhoff  result  for 
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u*  which  is  deduced  from  the  result  for  the  Kirchhoff  diffraction  by  a 
2  k 

conducting  half-plane.  The  u-|  term  is  now  viewed  as  being  the  h.f.  field 
diffracted  from  the  edge  P,  of  a  ha If -plane  located  at  x=0  for  -*<y^a. 
Similarly,  the  u£  term  is  viewed  as  being  the  diffracted  field  from  the 
edge  of  a  half  plane  located  at  x=0  for  -a<y<*.  In  doing  so,  one  assumes 
that  the  h.f.  fields  diffracted  by  the  strip  in  Fig.  5a  is  the  superposition 
of  the  fields  associated  with  rays  diffracted  from  the  edges  P^  and  P^  of 
two  overlapping  half  planes  (with  their  common  domain  of  overlap  =  width 

of  the  strip);  this  assumption  is  reasonable  since  ka  is  assumed  large.  The 

Kirchhoff  result  for  the  field  diffracted  by  the  edge  P^  of  a  half-plane 

when  a  plane  wave  is  normally  incident  on  it  is  given  in  terms  of  a  Kirchhoff 

edge  diffraction  coefficient,  as 


u{(P)  *  u1  (P1 )  DkU) 


-ik(p-asin$) 


where  u1^)  is  unity,  and 

(23)  Dk(*)  »  -  \  eW4  (-  FIX) 

with 

00 

2 

(24a)  F(X)  s  2i  l^e' lxI  J  e''T  dx  • 

i-TI 

and 

(24b)  X  =  2kpcos2  £  . 


L 

u^P)  of  Eq.  (22)  with  Dk (4> )  as  in  Eq.  (23)  approaches  the  following 
result  as  <p+rr: 


F  (X )  , 
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which  is  identical  to  the  first  term  on  the  RHS  of  Eq.  (21)  except  for 
the  additional  factor  F(X). 

F ( X )  %  (4x  -  2Xei7r/4)  •  ei(rr/4  +  X) 

as  <fmr ,  thereby  ensuring  a  bounded  result  for  u-j(P)  of  Eq.  (22)  at  the  shadow 

boundary  correspondi ng  to  The  derivation  of  the  result  in  Eq.  (22) 

is  summarized  in  Appendix  II.  The  function  F(X)  was  introduced  by 

Hutchins  and  Kouyoumji an[10]  to  describe  the  h.f.  diffracted  fields  in 

the  transition  regions  of  the  shadow  and  reflection  boundaries  associated 

with  a  perfectly  conducting  wedge  illuminated  by  a  plane  wave.  A  plot 

k 

indicating  the  behavior  of  F(X)  is  given  in  Fig.  6.  The  field  u^P)  is 
similarly  expressed  in  terms  of  as 

k  a  -i  k  (p+asi  n<j> ) 

(26)  u*(P)  *  u'(P2)  Dk(27T-4»)  ^ - 

/p 

which  is  also  finite  at  $=ir.  The  exponential  phase  factors  e±1'kasin,l>  -jn 
Eqs.  (22)  and  (26)  arise  from  choosing  the  origin  as  the  phase  reference. 

The  results  in  Eqs.  (22)  and  (26)  facilitate  the  calculation  of  the 

1/  U 

Kirchhoff  diffracted  fields  u-j  and  u^  over  a  greater  angular  region  about 
the  shadow  boundary  than  just  within  the  immediate  neighborhood  of  the 
shadow  boundary  at  <t>=ir.  It  can  be  shown  by  a  careful  limiting  operation 
that  a  superposi tion  of  Eqs.  (22)  and  (26)  also  leads  to  the  result  on  the 
LHS  of  Eq.  (21)  for  ^tt.  In  this  limiting  operation,  the  fields  in  Eq.  (22) 
and  Eq.  (26)  are  initially  assumed  to  emanate  from  P^  and  P^,  respectively 
to  the  far-zone  field  point  P  (in  the  vicinity  of  4>=it)  along  nearly  parallel 
ray  paths  which  converge  at  P.  Next,  the  field  corresponding  to  this  nearly 
parallel  convergent  ray  system  is  allowed  to  approach  the  parallel  ray 
system  (of  Fig.  5a)  in  the  mathematical  limit.  One  may  thus  employ  the 
field  representations  of  Eqs.  (22)  and  (26)  for  the  u-|  and  U£  terms  present 
in  Eqs.  (19)  and  (20),  respectively,  to  calculate  the  Kirchhoff  diffraction 

contribution  to  the  scattered  field  u|  over  a  larger  jngular  domain  about 

h  s 

4>=7t .  However,  before  the  results  of  Eqs.  (19)  and  (20)  for  us  can  be 

h 
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extended  to  calculate  the  scattered  field  over  this  larger  angular  domain 
about  4»=ir ,  one  must  modify  the  terms  involving  p*U)  and  q*(0  (which  are 
viewed  as  corrections  to  Ki rchhoff  diffraction)  in  a  manner  corresponding 
to  that  employed  for  U]  in  Eqs.  (22)  and  (26).  From  Eqs .  (A -1 0 )  and  (A-24) 
of  Appendix  I:  2 


(27a)  ?U)  =  L  jf  .-!«•  ♦ 


where 


(27b) 


.  p*(5)  if  the  operator  Q  =  1 

F(«)  =  i  3 

q*(0  if  the  operator  Q  =  ^ 


and  5  =  4-j  or  The  function  QV(t)  is  directly  related  to  the  Miller 
type  Airy  function  as  described  in  Appendix  I.  Clearly,  the  K  rchhoff 
correction  terms  in  Eqs.  (19)  and  (20)  may  be  expressed  in  terms  of 
F^U)  in  Eqs.  (27a)  and  (27b)  as: 

(28a)  -  e-'^4  (j£)i/3  ?(«)  e'11'4*  = 

-  -  *-”/4  (F 4'ika*  ^  Jr;-UT  - 

-  (unbounded  Ki rchhoff  result  of  Eq.  (21)  for  $%it), 

with  5  =  C]  corresponding  to  =  i|»i ,  as  before.  The  unbounded  Ki  rchhoff 
result  for2ij>%ir  is  seen  from  Eq.  (2?)  to  be: 


(28b) 


sika($-w) 


•  - 


The  integral  over  the  contour  lj  in  Eq.  (27a)  is  unbounded  at  4=0  (which 
corresponds  to  the  shadow  boundary  at  ♦=»),  but  f^(t)  in  Eqs.  (27a)  and 
(27b)  (and  therefore  the  result  in  Eq.  (28a))  is  bounded[4,7]  because  the 


19 


singularity  in  the  — —  term  of  Eq.  (27a)  exactly  cancels  the  singularity 
2/^5 

in  the  integral  (over  Ip  when  £->0.  It  is  then  obvious  that  the  Kirchhoff 
correction  terms  in  Eqs.  (19)  and  (20)  can  be  appropriately  modified  to  be 
useful  over  a  larger  angular  domain  about  the  shadow  boundary  (<|>=ir),  if  the 
unbounded  Kirchhoff  result  for  4>  in  the  close  vicinity  of  ir  in  Eq.  (28a)  is 
replaced  by  the  unbounded  Kirchhoff  result  valid  over  a  larger  angular 
domain  about  This  is  accomplished  by  noting  from  Eq.  (22)  that  the 

factor  F(X)  keeps  Dk(<}>)  bounded  at  $=ir,  and: 


bounded  Kirchhoff 
diffraction  field 


unbounded  Kirchhoff 
diffraction  field 


x  F(X) . 


Thus,  Eq.  (28a)  which  is  valid  for  is  now  modified  by  replacing 

-ikai|>i 

-e-W4  |T  (|1),/3?U))e 

2 


(28c)  -  e"W4  (ibjr)173 


-ikaip] 


2  e-ikp 


?Ui) 

2 


1 

2*4^1 


^PT  „-1,/4  ±1  e-iK(«asin») 

i2ir¥  e  si  n4>  r 


Combining  the  Kirchhoff  diffraction  results  of  Eq.  (22)  and  Eq.  (26) 
with  the  Kirchhoff  correction  terms  of  Eq.  (28c)  yield  an  expression  for 
u|(P)  which  is  useful  over  a  larger  angular  domain  about  $=*  (in  contrast 

to  Eq.  (19)  and  Eq.  (20)  which  are  good  only  in  the  close  vicinity  of  $=ir): 


(29a) 


uS(P,P1)  *  - 
2 


e-iw/4 

±sinp 


[F(X] )-l ] 
2 


J-ik(p+asinp) 

/p  ) 


( 


1/3 


V 


-ika^i 


1 

2  e-W4 


1 


L 


f<V  ■  ii 


-iks 


*S 


with  X-j  =  2kp  cos24>/2,  and  =  2kp  cos2  (2ir-<j>  )/2 . 

uS (P ,P-j )  corresponds  to  the  scattered  field  at  P  arising  from  the 
2 

penumbral  region  associated  with  or  P^.  Also 


u^(P,P,),  when  =  P*(ti),  for  the  TM 


(29b)  uS(P,P1 )  =  ^ 

2 


's'’  ’■  1 
2 


uj^P.Pj),  when  ^(^)  =  q*(€1),  for  the  TE 


The  total  h.f.  field  at  a  point  P  arising  from  the  vicinity  of  P^  on  the 
cylinder  may  be  written  as: 


(30a) 


us(P,P1)  * 


u’(P>  +  u'(P,P,) 
USS(P,P,) 


for  the  TMz  case , 


on  the  lit  side  (4><ir) 
on  the  shadow  side  U>n) 


case 

case. 


and 


r 


(30b)  u^P.P,)  ^ 


uj(P)  +  u^P.Pp 

uh<P’V 


for  the  TEz  case. 


on  the  lit  side  (4><w) 
on  the  shadow  side  {$>■*) 


The  total  h.f.  field  at  P  arising  from  the  vicinity  of  on  the  cylinder 
is  similarly  expressed  as: 


21 


on  the  lit  side  (4><w ) 
on  the  shadow  side  (<(>>Tt) 


on  the  lit  side  (<j><ir) 
on  the  shadow  side  ) . 

for  the  T£z  case. 

The  total  h.f.  field  at  P  in  the  transition  region  is  the  sum  of  the 
fields  arising  from  the  vicinity  of  and  P 2  on  the  cylinder,  thus: 

(31a)  u$(P)  =  us(P,P1)  +  us(P,P2)  for  the  TMz  case 

and 

(31b)  uh(P)  =  uh(P,P1)  +  uh(P,P2)  for  the  TEz  case. 

The  results  of  Eqs.  (31a)  and  (31b)  will  be  employed  in  part  B  which 
follows,  to  calculate  the  fields  in  the  transition  region  associated 
with  the  canonical  problem  of  Fig.  2.  One  notes  that  the  results  in 
Eqs.  (29a)  and  (29b)  are  based  on  a  direct  generalization  of  the  results 
in  Eqs.  (19)  and  (20).  The  accuracy  of  these  results  will  be  indicated  in 
the  next  section  wherein  they  are  applied  to  specific  problems.  Possibly, 

a  more  accurate  asymptotic  h.f.  result  for  us(P,P-|)  (when  P  is  not 

2 

necessarily  in  the  close  vicinity  of  shadow  boundary)  obtained  from 
first  principles  might  involve  terms  that  are  not  a  simple  superposition 
of  the  effective  edge  diffracted  field  term,  and  a  surface  dependent  field 
term  as  in  Eq.  (29a);  such  an  asymptotic  solution  based  on  mathematical 
rigor  rather  than  on  simplifying  physical  assumptions  to  bypass  mathematical 
complexities  appears  to  be  a  difficult  task.  However,  judging  from  the 


(30c)  us(P,P2)  <v 


and 


(30d) 


u’ (P)  +  U^(P,P2) 
Uf(P.P2) 


for  the  TMZ  case. 


uh(P.P2)  * 


uj(P)  +  ujj(P,P2) 
u^(P,P2) 
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numerical  results  obtained  from  the  use  of  Eq.  (29a)  in  the  following 
sections,  it  appears  that  the  present  solution  yields  sufficiently  accurate 
values  for  engineering  applications. 


B.  H.F.  Scattering  by  a  Perfectly-Conducting 
2-b,  Rounded-End  Thick-Screen 


The  canonical  configuration  of  interest  is  illustrated  in  Fig.  2. 

It  is  desired  to  calculate  the  total  h.f.  field  surrounding  this  2-D 
rounded-end,  thick-screen  configuration  when  it  is  excited  by  a  normally 
incident  electromagnetic  plane  wave  of  the  TMz  or  the  TEz  type.  The 
solution  to  this  canonical  problem  is  of  relevance  to  a  study  of  the 
effects  of  the  aircraft  fuselage  on  the  radiation  from  wing  mounted  antennas, 
as  indicated  previously.  Figure  2  illustrates  the  junctions  and 
resulting  from  the  smooth  join  of  the  half-cylinder  onto  the  thick  half 
plane.  In  the  present  analysis,  the  diffraction  effects  arising  from  the 
junctions  J-j  and  J2  are  not  included  for  the  following  reasons: 

1)  The  actual  fuselage  of  an  aircraft  will  be  approximated  by 
this  model  for  which  these  junctions  will  only  be  artificially 
added  to  our  model  in  that  they  are  not  present  in  the  actual 
aircraft  structure. 

2)  The  junction  effects  are  important  only  if  the  penumbral 
region  lies  in  the  neighborhood  of  the  junction,  or  if  the  field 
point  lies  in  the  neighborhood  of  the  reflection  boundary  associated 
with  the  junction.  It  is  clear  therefore,  that  in  the  present 
problem,  which  deals  with  normal  incidence  and  with  aspects  within 
and  near  the  transition  region,  the  junction  effects  are  negligible. 

The  total  h.f,  field  in  the  lit  region  for  the  problem  in  Fig.  2  is 
adequately  described  by  the  fields  associated  with  the  incident  and 
reflected  rays  of  geometrical  optics,  whereas,  the  h.f.  field  in  the  shadow 
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region  is  accurately  described  by  the  surface  diffracted  rays  of  Keller's 
geometrical  theory  of  di ffraction[9]  (hence  by  the  GTD).  However,  the 
geometrical  optics  field  description  fails  near  and  on  the  reflection 
boundary  (which  is  also  the  shadow  boundary  in  this  case).  On  the  other 
hand  the  GTD  field  description  requires  an  increasing  number  of  surface  ray 
modes  to  accurately  describe  the  field  near  the  shadow  boundary  thereby 
making  the  GTD  result  cumbersome  to  employ.  Thus,  within  and  near  the 
transition  region  adjacent  to  the  shadow  boundary,  one  may  employ  the 
result  developed  in  part  A  to  calculate  the  total  h.f.  fields  there.  The 
specific  results  for  the  lit,  transition,  and  shadow  regions  associated 
with  the  canonical  configuration  of  Fig.  2  are  indicated  below. 

Field  Description  in  the  Lit  Region  -  The  total  field  in  the  lit 
region  is  a  superposition  of  the  incident  field  u1 ,  and  the  reflected  field 
ur  predicted  by  geometrical  optics,  in  which 


(32a) 


u(  *  4  -  e1k\ 
h 


for  TMz  case. 


(32b) 


ur  =  u£  =  u^(PD)  •  R 


S1|pr+s 


where  the  subscript  s  is  used  for  the  TMz  case,  and  the  subscript  h  is 
used  for  the  TEz  case  as  before.  The  result  in  Eq.  (32b)  is  in  a  ray- 
optical  form[ll]  which  is  valid  in  the  near  zone.  PR  denotes  the  point 

of  reflection  in  Fig.  7a,  and  Uj(PR)  is  the  value  of  the  incident  (jez) 

h  Z 

field  at  PR.  The  reflection  coefficient,  Rs  is  defined  by: 

h 


(32c) 

(32d) 


Rs  =  -1, 


Rh  =  +  1  * 


for  the  TMz  case. 


for  the  TEZ  case. 


r 


Fig.  7a.  Reflected  ray  geometry. 

s  denotes  the  distance  from  PR  to  the  field  point  P  in  the  lit  zone.  The 
quantity  within  the  radical  in  Eq.  (32b)  refers  to  the  caustic  distance 
(in  the  plane  of  reflection)  associated  with  the  reflected  wavefront;  pr 
is  found  to  be: 

a  cos* 

(32e)  Pr  =  - 2“-^  . 


The  total  field  at  a  point  P  in  the  lit  zone  is  thus  given  by: 

-iks 


(32f)  us(P)  *u£(P)  +  uJ(PR)  (*1) 
h  h  h 


a  cos*. 


TM_ 


a  cos*Q  +  2s 


;  for  case. 
Lz 
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If  far  zone  results  are  desired,  one  must  enforce  the  condition  s  »  p 

r 

in  Eq.  (32b).  The  result  in  Eq.  (32f)  is  valid  in  the  lit  zone  for  aspects 
sufficiently  far  from  the  reflection  boundary  (which  occurs  when  P^  coincides 
with  P^ ). 

Field  Description  in  the  Transition  Region  -  The  total  field  in  the 
transition  region  may  be  immediately  written  down  from  Eqs.  (30a)  and  (30b) 
given  in  Part  A.  It  is  conjectured  here  that  the  fields  in  the  transition 
region  of  Fig.  2  are  identical  to  those  in  Eqs.  (30a)  and  (30b)  which 
arise  from  the  penumbral  region  associated  with  P-j  on  the  circular  cylinder 
of  Fig.  1.  This  conjecture  is  based  on  the  h.f.  approximation  which  allows 
one  to  consider  h.f.  diffraction  as  being  a  local  phenomenon.  The 
results  in  Eqs.  (30a)  and  (30b)  for  the  TMz  and  the  TEz  cases,  respectively, 
are  valid  in  the  far  zone;  however,  since  h.f.  diffraction  is  a  local 
phenomenon,  one  may  also  directly  modify  the  far-zone  results  of  Eqs. 

(30a)  and  (30b)  to  obtain  the  following  near  zone  transition  fields  at  P: 


with  the  understanding  that  f^(0  =  P*(0  for  u$(P)  and  fye)  =  q*U)  for 
u^(P);  also  H(tt-o)  is  the  Heaviside  step  function  whose  value  is  unity  if 
a<ir  and  zero  if  a>ir.  The  near  zone  result  of  Eq.  (33)  is  obtained  by 
simply  replacing  the  far  zone  ray  distances  by  their  corresponding  near 
zone  ray  distances.  The  quantities  s,  $  and  5  in  Eq.  (33)  are  given  in 
terms  of  the  near  field  coordinates  by: 


(34a) 


s 


in  the  lit  portion  of  the  transition  region 
in  the  shadowed  portion  of  the  transition  region 


i  -  i 4* £ | »  in  the  lit  portion  of 
(34b), (34c)  5  =  (^-)1/3  4,  ;  *  =  i  the  transition  region. 

'  j<l>  j ,  in  the  shadowed  part  of 

the  transition  region. 

where  the  angles  <i> ^  and  ips,  and  the  distances  sjL  and  ss  are  indicated  in 
Fig.  7b  below.  The  F(X)  function  appearing  in  Eg.  (33)  is  defined  earlier 


P  (  LIT) 


Fig.  7b.  Equivalent  ray  geometries  for  the  fields  in  the  lit  and 
shadowed  portions  of  the  transition  region  associated 
with  the  shadow  boundary  at  P-j . 


in  Eq.  (24a).  The  new  argument  (for  near  zone  calculations)  of  F(X) 
is  clearly: 


(35a) 


X  =  2k d  cos2  S- 


where  a  is  the  angle  (measured  in  the  cw  sense)  subtended  by  r.  Kirchhoff 
diffracted  ray  (directed  along  P^P)  and  the  negative  of  the  incident  ray 
direction  (i,e.,  along  +  x). 


(35b) 


in  the  lit  part 
in  the  shadowed  part. 


The  distance  d  is  measured  from  P^  to  P  and 


(35c) 


d  = 


in  the  lit  part 
in  the  shadowed  part. 


The  quantities  a^,  ,  d^  and  dg  are  illustrated  in  Fig.  7b. 


The  transition  region  solution  discussed  above  must  serve  a  two¬ 
fold  objective,  namely  that  it  should  provide  a  sufficiently  accurate  h.f. 
field  description  within  the  transition  region;  secondly,  it  should  blend 
smoothly  with  the  geometrical  optics  solution  in  the  lit  region,  and  also 
blend  smoothly  with  the  GTD  solution  in  the  shadow  region.  Some  numerical 
calculations  which  indicate  the  accuracy  of  the  present  results  are  provided 
later. 


Field  Description  in  the  Shadow  Region  -  For  aspects  in  the  deep 
shadow  region,  the  scattered  field  (the  incident  field  is  zero  within  the 
shadow  region)  is  given  a  ray  optical  representation  valid  in  the  near 
zone.  The  result  based  on  GTD  is  given  as[ll]: 
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since  (P-j )  =  1.  The  geometrical  details  are  given  in  Fig.  7c  for  the 
h 
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Fig.  7c.  Ray  geometry  for  the  GTD  fields. 

surface  ray  modes  which  traverse  a  path  P-jPpP  to  the  field  point.  Even 

though  the  inclusion  of  an  infinite  number  of  GTD  surface  ray  modes  is 

indicated  in  Eq.  (36),  it  is  generally  found  that  sufficiently  accurate 

results  for  the  field  are  obtained  by  summing  only  the  first  five  modes 

(corresponding  to  p=l  to  p=5),  provided  the  field  point  P  is  far-removed 

s  s 

L  L 

from  the  transition  region.  Dn  and  are  the  cylinder  diffraction  coef- 

P  P  .  . 

ficient  and  attenuation  constant  for  the  p-th  surface  ray  mode,  respectively. 

D.R.  Voltmer  (in  his  Ph.D  dissertation,  Ohio  State  University,  1970)  has 

obtained  diffraction  coefficients  and  attenuation  constants  to  order 

(|-^—  )^3  (here  p  =  local  radius  of  curvature  of  the  surface  along  the  ray 

path)  which  are  valid  for  smooth  convex  surfaces  of  non-constant  curvature 
(this  includes  the  circular  cylinder,  and  the  sphere  as  special  cases);  a 


table  of  his  diffraction  coefficients  and  attenuation  constants  is,  also, 
available  in  Reference  [8].  The  superscripts  s  and  h  on  the  diffraction 
coefficient  and  attenuation  constant  denote  the  TMz  and  the  TEz  cases, 
respectively. 

Numerical  Results  -  The  numerical  results  for  the  near  field  patterns 
of  a  circular  cylindrically  tipped  thick  half  plane  (as  shown  in  Fig.  2) 
illuminated  by  a  TM_,  (acoustic  soft)  or  TEz  (acoustic  hard)  plane  wave  have 
been  calculated  via  the  formulas  developed  above.  These  results  are 
illustrated  in  Figs.  8  through  10.  The  geometrical  optics  field  shown 
as  a  dot-dashed  line,  consists  of  the  incident  and  reflected  fields  up  to 
the  shadow  boundary.  The  transition  field,  shown  as  a  solid  line,  starts 
at  <t>=90°  and  is  shown  into  the  deep  shadow  region.  The  surface  ray  mode 
field  (obtained  by  summing  over  only  the  first  5  modes)  is  shown  as  a 
dotted  line  and  is  valid  in  the  deep  shadow  region.  Note  that  the 
geometrical  optics  field  and  the  transition  field  blend  well  in  the  angular 
range  90°<<J><1 20° .  The  transition  field  and  the  surface  ray  mode  field  for 
the  hard  cases  blend  in  the  neighborhood  of  $=180°.  The  soft  cases  do  not 
blend  as  quickly.  The  phase  for  the  above  near  field  patterns  also  blended 
well;  however,  it  is  difficult  to  illustrate  and  will  not  be  shown  here. 

The  validity  of  the  above  high  frequency  solutions  is  demonstrated 
by  applying  them  to  the  problem  of  plane  wave  scattering  by  a  circular 
cylinder  (as  shown  in  Fig.  1).  The  results  are  then  compared  with  the 
exact  solution  for  the  problem  which  is  found  from  the  eigenfunction 
expansion  in  Figs.  11a  and  lib.  The  high  frequency  solution  is  composed 
of  the  superposition  of  fields  from  the  two  penumbral  regions  associated 
with  P-|  and  ?£  using  the  appropriate  results  (of  part  B)  in  their 
respective  regions  of  validity.  The  good  agreement  lends  confidence  to 
the  present  high  frequency  solutions. 
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Fig.  8a.  Near  field  pattern  for  cylindrically  tipped 
thick  half-plane,  hard  case,  ka=10,  kp=100. 
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.  8b.  Near  field  pattern  for  cylindrically  tipped 
thick  half-plane,  soft  case,  ka=10,  kp=100. 
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Fig.  9a.  Near  field  pattern  for  cylindrical ly  tipped 
thick  half-plane,  hard  case,  ka=20,  kp=50. 
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Fig.  9b.  Near  field  pattern  for  cylindrically  tipped 
thick  half-plane,  soft  case,  ka=20,  kp=50. 
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Near  field  pattern  for  cylindrically  tipped 
thick  half-plane,  hard  case,  ka=40,  kp=100. 
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Fig.  11a.  Comparison  of  high  frequency  solution  with 
exact  solution  for  near  field  pattern  of 
circular  cylinder,  hard  case,  ka=20,  kp=50. 


11b.  Comparison  of  high  frequency  solution  with 
exact  solution  for  near  field  pattern  of 
circular  cylinder,  soft  case,  ka=20,  kp=50. 


III.  HIGH  FREQUENCY  SCATTERING  BY  A  3-D, 

SMOOTH  CONVEX  SURFACE 

As  developed  in  the  previous  section  the  transition  solution  can  be 
divided  into  two  terms.  These  are  the  Kirchhoff  term,  which  is  independent 
of  the  polarization  and  curved  surface  properties,  and  a  second  term  which 
is  dependent  on  the  surface  curvature  and  associated  electrical  properties. 
The  usefulness  of  this  solution  is  that  each  of  these  terms  can  be  system¬ 
atically  modified  to  three-dimensional  surfaces  (based  on  the  high  frequency 
assumption  that  diffraction  is  a  local  phenomenon).  One  notes  that  the 
geometrical  optics  and  geometrical  theory  of  diffraction  solutions  are 
well-known  for  the  three-dimensional  case  as  presented  in  Reference  [11]. 

In  this  section,  the  near  zone  scattering  by  a  spherically  capped 
semi-infinite  circular  cylinder  when  excited  by  a  normally  incident  plane 
wave  is  examined  as  shown  in  Fig.  12.  The  object  of  this  study  is  to 
illustrate  the  three-dimensional  nature  of  our  high  frequency  solution  and 
to  verify  that  it  blends  nicely  with  the  geometrical  optics  field.  The 
field  surrounding  the  configuration  in  Fig.  12  is  given  below  in  the  lit, 
transition,  and  shadow  regions  in  a  manner  similar  to  that  indicated  for 
the  canonical  configuration  of  Fig.  2,  in  part  B  of  the  previous  section. 
Once  again,  the  diffraction  effects  arising  from  the  junction  of  the  hemi¬ 
spherical  cap  and  the  cylinder  are  ignored  for  reasons  indicated  previously. 


Field  Description  in  the  Lit  Region  -  The  geometrical  optics  solution 
composed  of  the  incident  field  (u1 )  and  the  reflected  field  (ur)  observed 
at  a  near  zone  field  point  P  in  the  lit  zone  is  explicitly  given  by: 


The  point  of  reflection  PR  and  other  geometrical  parameters  (such  as  a,  s, 

♦0,  etc.)  are  the  same  as  in  Fig.  7a.  The  reflection  coefficient  Rs  was 

h 


39 


P<  LIT) 


Fig.  12.  Geometrical  configuration  for  the  hemi spherically 
capped  semi -infinite  circular  cylinder. 

previously  defined  in  Eqs,  (32c)  and  (32d).  Ug(PR)  denotes  tfie  value  of 

ul  at  the  point  of  reflection  PD  on  the  hemispherical  cap. 
h  R 

Field  Description  in  the  Transition  Region  -  The  total  field  in  the 
transition  region  for  the  3-D  geometry  of  Fig.  12  can  be  obtained  by 
generalizing  the  2-D  results  for  the  geometry  in  Fig.  2  which  are  indi¬ 
cated  in  Eq.  (33).  This  generalization  from  2-D  to  3-D  is  based  on  the 
ray-optic  description  for  the  fields;  analogous  generalizations  have  been 
previously  indicated  elsewhere[8,9,ll  ,12]  for  certain  antenna  and  scattering 
problems.  We  obtain 
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(38a)  us(PK(P)H(,-a)+{-^e-1^  ^1H=±±  j e 


-({£  (j£)1/3  e-ika*  e"iir/4 


where  H(ir-a),  f^U)>  F(X),  a,  s,  d  and  have  the  same  meaning  as  in  Eq. 

(33).  The  quantity  f(s)  is  given  by: 
f 


\ 

fte)  -  ~~ 

e'ik$  f(s) 

_  24c__ 

! 

(38b) 


f(s)  =  /  ^s(p+s) 


in  the  illuminated  region  (a<n) 


dn  - 

Idn  i sTp+sJ  in  the  shadow  region  (a>ir). 


The  geometrical  details  are  the  same  as  those  indicated  in  Fig.  7b.  The 

dn 

quantity  is  the  spatial  spread  of  the  field  on  the  surface  of  the 

hemispherical  cap  along  the  ray  path  P^P$  of  Fig.  7b.  In  particular,  dnQ 

is  the  width  of  two  adjacent  rays  at  P, ;  whereas,  dn  is  the  width  of  the 

dn 

same  two  adjacent  rays  evaluated  at  P  .  The  quantity  ^gp  is  given  by: 


(38c) 


i  dnn 

i&r 


N  sine 


where  0  = 


0 1 i t  in  the  lit  penumbra  region  (=tt/2  -  |<|»A!)  i.e.,  for  e<| 
esh  in  the  shadow  penumbra  region  (=ir/2  +  ks  I )  i-e.,  for  e> 


Also,  the  caustic  distance  p  appearing  in  Eq.  (38a)  is  the  radius  of 
curvature  of  the  diffracted  wavefront  in  the  plane  _l  <-y  plane,  but  contain¬ 
ing  the  ray  shed  from  the  surface;  it  is  given  by 


(38d) 


p  =  atane. 


Incorporating  Eqs,  (38c)  and  (38d)  in^^p  ^-(p-Ppy  of  Eq.  (38b)  leads 


to 
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(38e) 


Thus,  p>0  when  8=6^ t>  whereas  p<0  for  e=esh-  If 
sTp+s)  '  °’ 

then 


I 


s(p+s) 


1 


P 

s(p+s) 


.iff/ 2 


because  a  caustic  has  been  crossed.  One  notes  that  bas  been  replaced 
by  unity  in  the  lit  part  of  the  transition  region  (see  Eq.  (38b))  because 
it  is  conjectured  that  the  fields  scattered  into  the  lit  portion  should 
be  independent  of  the  spatial  spread  factor  of  the  rays  along  on  the 
surface  (this  fact  is  directly  evident  in  the  reflected  ray  description  for 
the  lit  region). 


In  order  to  keep  the  Kirchhoff  diffracted  term  bounded  one  notes 
that  the  X  in  F(X)  of  Eq.  (38a)  should  be 


(38f ) 


X  =  2k 


d(pe+d) 


cos2a/2. 


The  quantity  pg  is  the  caustic  distance  for  the  Kirchhoff  diffracted  ray 
in  the  plane  containing  the  diffracted  ray  and  the  edge  of  the  curved 
screen  at  P^;  it  is  computed  easily  via  the  formula  for  caustic  distances 
for  curved  edges  given  by  Kouyoumjian[12]. 


Field  Description  in  the  Shadow  Region  -  The  <  iffracted  ray  field  with¬ 
in  the  shadow  region  at  a  field  poin'-  P  in  the  near  zone  is  obtained  by 
extending  the  2-D  result  of  Eq.  (36)  to  the  3-D  result  in  the  usual  manner[ll]. 
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(39)  us  (P)'MJ?  (P)  = 
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p=l  h  K 
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P=1  P 


Dp(fV  1s(asine+scos0)  e 
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where 


4(9 })  =  !> 

h 


and 


dn 

dn  \ 

p  -  i 

s(p+s) 

s (asi  n8+scose) 1 


as  noted  previously  in  Eq.  (38e).  The  diffraction  coefficient  D*1,  and  the 

s  P 

attenuation  constant  are  indicated  in  Reference  [8]  for  the  spherical 
surface.  The  geometrical  details  pertaining  to  the  GTD  ray  description  of 

Eq.  (39)  for  us(P)  are  the  same  as  those  in  Fig.  7c.  If 
h 

sTp+s)  °- 

fdnl 


then  IdTrlii^?) 

must  be  interpreted  as 
Idn, 


J*/2 


due  to  the  crossing  of  a  caustic. 
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It  is  evident  from  the  results  in  Eqs.  (38a),  (38b),  and  (39)  that 
a  caustic  of  the  diffracted  rays  is  present  for  4>=tt  (and  a=-sin"^  ^  in 
Eq.  (38a)).  One  notes  that  on  crossing  the  caustic  which  occurs  at  4>=Tr 
(i.e.,  along  the  negative  x-axis), two  additional  diffracted  rays  which 
creep  around  either  side  of  the  cylindrical  surface  also  contribute  to  the 
field  at  P;  the  field  of  these  rays  is  easily  calculated  by  GTD.  The 
caustic  field  may  be  assumed  to  result  from  the  radiation  by  equivalent 
half-ring  electric  and  magnetic  currents  in  which  the  half-ring  is  formed 
by  the  locus  of  the  points  of  diffraction  of  the  rays  shed  from  the 
hemispherical  cap  to  the  field  point  P  for  which  4>=ir.  The  equivalent  ring 
current  concept[8,14,15,16]  proves  to  be  a  useful  tool  for  evaluating 
the  proper  field  behavior  at  and  in  the  vicinity  of  a  caustic.*  Such  an 
alternative  field  representation  for  the  caustic  region  is  necessary 
because  the  simple  field  description  of  Eqs.  (38ab)  and  (39)  fails  in 
this  region  (the  results  of  Eqs.  (38a)  and  (39)  become  infinite  at  the 
caustic).  The  transition  from  the  single  diffracted  ray  arriving  at  P 
in  the  shadow  region  (whose  field  is  given  by  Eq.  (39))  before  crossing  the 
caustic  to  the  three  diffracted  rays  which  arrive  at  P  in  the  shadow  region 
after  crossing  the  caustic  (the  three  rays  result  from  the  ray  corresponding 
to  Eq.  (39)  shed  from  the  hemispherical  cap  to  P  plus  the  two  additional 
rays  which  creep  around  either  side  of  the  cylindrical  surface  before 
being  shed  towards  P)  will  not  be  treated  in  the  present  report.  This 
transition  is  described  by  the  solution  obtained  from  the  equivalent 
current  concept  for  the  caustic  region.  The  caustic  field  analysis  will 
be  completed  in  a  future  report. 

Numerical  Results  -  The  near  field  patterns  for  the  hemispherically 
capped  cylinder  of  Fig.  12  for  a  radius  of  ka=10  and  observation  distance 
of  kp=100  are  presented  in  Figs.  13a  and  13b,  when  this  structure  is 
illuminated  by  a  TEz  or  TMz  type  plane  wave,  respecti vely.  The  various 
fields  are  represented  as  done  previously  for  the  2-D  problem.  As  can  be 
seen  the  transition  solution  blends  well  with  the  geometrical  optics  field 
for  both  the  hard  and  soft  cases.  The  presence  of  the  caustic  at  4>=tt  , 

*"A"  study  of  the  caustic  due  to  the  curved  screen  Kirchhoff  diffraction  term 
(i.e.,  for  a=sin-‘  a/d  in  Eq.  (38a)  when  4>=w)  will  also  be  completed  in  a 
future  report. 
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Fig.  13a.  Near  field  pattern  of  hemispherically  capped 
cylinder,  hard  case,  ka=10,  kp=100. 
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however,  interferes  with  the  blending  of  the  transition  solution  with  the 
field  of  the  surface  ray  modes.  As  stated  above,  the  singularity  due  to 
the  caustic  can  be  easily  corrected;  however,  this  has  not  been  done  here 
but  will  be  completed  in  a  future  report.  When  this  caustic  correction  is 
incorporated,  it  is  expected  that  the  blending  of  the  transition  solution 
with  the  surface  ray  GTD  solution  will  be  as  good  as  that  indicated  earlier 
for  the  2-D  problem. 

IV.  CONCLUSIONS 

The  object  of  this  report  has  been  to  develop  the  complete  solutions 
for  the  fields  scattered  by  curved  surfaces  illuminated  by  near  zone 
antennas.  Specifically  the  high  frequency  solutions  for  curved  surfaces 
have  been  completed  for  several  years  except  for  a  transition  solution 
which  is  needed  around  the  shadow  boundaries.  An  approximate  transition 
region  asymptotic  solution  is  presented  and  developed  in  such  a  form  that 
it  can  be  extended  to  arbitrary  curved  surfaces.  Combining  this  solution 
with  the  geometrical  optics  term  in  the  lit  region,  and  the  geometrical 
theory  of  diffraction  term  in  the  shadow  region,  should  provide  a  complete 
high  frequency  solution  for  the  scattered  fields  of  antennas  illuminating 
curved  surfaces  in  the  near  zone. 

The  solution  is  presented  in  terms  of  the  thick  half  plane  with 
a  cylindrical  cap  showing  how  the  various  terms  blend  together  in  each 
region.  The  transition  region  solution  developed  here  (Eqs.  (33)  and 
(38a))  can  be  shown  to  analytically  reduce  to  the  GTD  surface  ray  mode 
solution  as  F(X)-*T  in  the  shadow  region  sufficiently  far  from  the  shadow 
boundary.  We  have  not  yet  shown  that  this  transition  solution  can 
analytically  reduce  to  the  geometrical  optics  reflected  field  as  F(X)-*1 
in  the  illuminated  region  sufficiently  far  from  the  shadow  boundary. 
However,  the  numerical  results  presented  indicate  that  the  transition 
region  solution  does  indeed  blend  very  nicely  with  the  GTD  and  the 


geometrical  optics  solutions.  Although  the  present  transition  region 
solution  appears  to  be  fairly  accurate,  it  would  be  desirable  to  obtain 
a  transition  region  solution  from  a  more  rigorous  asymptotic  analysis  of 
the  canonical  circular  cylinder  scattering  problem. 

Even  though  the  solution  presented  describes  the  one  term  as  being 
a  Kirchhoff  solution,  one  can,  also,  show  that  this  term  can  result 
using  the  GTD  edge  diffraction  coefficient  instead  of  the  Kirchhoff  edge 
diffraction  coefficient.  Using  this  approach  the  values  of  the  transition 
solution  do  not  change  much  near  the  shadow  region,  but  they  might  provide 
improvements  away  from  the  shadow  boundary.  This  approach  will  be  investi¬ 
gated  for  the  spherically  capped  cylinder  in  the  caustic  region  in  that 
in  this  sector  the  Kirchhoff  and  GTD  solutions  are  observably  different. 

While  the  transition  solution  is  being  completed,  the  solution  for 
a  finite  elliptic  cylinder  will  be  developed.  This  solution  will  provide 
a  versatile  model  for  the  fuselage  when  the  antenna  is  not  mounted  directly 
on  it,  i.e.,  such  as  wing  or  pod  mounted  antennas.  It  can,  also,  be  used 
to  examine  the  scattering  effects  of  stores,  etc.  With  the  completion  of 
this  solution,  one  should  have  all  the  necessary  theoretical  tools  to 
analyze  any  high  frequency  on -aircraft  antenna  pattern  using  a  reasonable 
model  representation  of  the  actual  aircraft. 


APPENDIX  I 


EVALUATION  OF  THE  INTEGRALS  FOR  THE 
FIELD  IN  THE  TRANSITION  REGION 


A  summary  of  the  details  concerning  the  asymptotic  evaluation  of 
the  integrals  in  Eqs .  (7a)  and  (7b)  of  Section  II  is  provided  in  this 
appendix.  The  TMz  case  (correspondi nq  to  Eq.  (7a))  is  considered  first; 
it  is  followed  by  the  TEz  case  (of  Eq.  (7b)). 


TM  case:  Employing  Eqs.  (9a)  and  (13)  in  Eq.  (7a)  leads  to: 
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then  Eq.  (A-l)  may  be  written  as: 
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with  the  understanding  that  only  the  *,=0  term  is  of  interest  for  the  reason 
mentioned  in  Section  II.  The  integral  Ij,  of  Eq.  (A-2)  is  evaluated  first. 
One  begins  by  incorporating  Eq.  (10a)  in  Eq.  (A-2);  thus: 
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with  =  (!p-)^3  .  The  sum  on  l  can  be  incorporated  in  the 

integration  if  one  changes  the  limits  of  integration  to  -•»<({> '<«>,  since  the 
stationary  points  occur  at  <t>  '=4>-(4A-l  )tt/2h4)^ .  Expanding  the  exponential 
e-ika[U  -u/2 )-cos )]  -n  ^  (a_5)s  about  the  stationary  point  <f>^ 

corresponding  to  i=0  only  (since  multiple  encirclements  are  neglected),  and 
retaining  terms  only  to  the  third  power  in  one  obtains: 
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after  re-arranging  terms,  and  interchanging  the  orders  of  integration. 
Introducing  a  change  of  variables: 

(A-7)  (jj2-)1/3(*'-+')  =  y  , 

allows  one  to  simplify  Eq.  (A-6)  to: 
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One  notes  that  the  integral  on  y  yields  the  Miller  type  Airy  function[13] 
denoted  by  Ai(-r);  specifically: 


i  r*  -i  (ty+  y 
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Defining  the  function  p*(s)  introduced  in  Secti  n  II,  namely[4]: 
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one  may  now  rewrite  Eq.  (A-8)  in  terms  of  Eq.  ^A-9),  and  Eq.  (A-10)  as: 
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where 

(A-12)  ^  =  (^-)1/3(4»-tt) 


A  similar  analysis  for  I2  in  Eq.  (A-3)  leads  to: 
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Incorporating  Eqs.  (A-l 1 )  and  ( A-l 3 )  in  Eq.  (A-4)  leads  finally  to: 
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where  and  *2  are  defined  in  Eqs.  (17a)  and  (17b).  The  result  in 
Eq.  (A-15)  corresponds  to  that  indicated  in  Eq.  (15). 

TE  case:  Employing  Eq.  (9b)  and  Eq.  (14)  in  Eq.  (7b)  leads  to: 
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then  Eq.  (A-16)  may  be  written  as: 
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Following  the  procedure  for  the  TMz  case,  one  may  write  ^  in  terms  of 
the  stationary  point  at  $g=*-w/2  (corresponding  to  the  4=0  term)  as: 


y  tm  -ika  (^-irJ-ijr— ($'-$' ) 

(A-20)  ^  %  J  g(€1)(#,-^)e  6  s 


d4> 1  , 


where  cos ($-$')  in  Eq.  (A-17)  has  been  approximated  by  ($'-4^)  in  the 
neighborhood  of  the  stationary  point,  and  the  limits  of  integration  have  been 
changed  to  !♦'!<"  as  done  for  the  TMZ  case.  Incorporating  Eq.  (10b)  into 
Eq.  (A-20)  gives: 
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Introducing  Eq.  (A-7)  into  Eq.  (A-21)  yields: 
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One  recognizes  the  integral  on  y  to  be  related  to  the  derivative  of  the 
Airy  function  Ai(r);  namely[13]: 


-i  (ty+  £  )  _ 

(A-23)  —  dYve  3  =  2ik7  Ai'(x)  ee  2i  V'(t). 

nr  -• 


One  also  defines  the  function  q*{£)  introduced  in  Section  II;  namely[4]: 
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Incorporating  Eqs.  (A-23)  and  (A-24)  in  Eq.  (A-22)  finally  gives: 
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with  £.|  defined  in  Eq.  (A-12).  A  similar  analysis  for  ^  leads  to: 


(A-26)  L  %  eika(*_1T^2i  (^-)’2/3  q*U9) - 1 


I  '• 


► 

f 


with  $2  defined  in  Eq.  (A-14).  Employing  Eqs.  (A-25)  and  (A-26)  in  Eq. 
(A-19)  gives 
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with  ij/-|  defined  in  Eq.  (17a)  and  Eq.  (17b).  The  result  in 
2 

corresponds  to  the  desired  result  indicated  in  Eq.  (16). 


Eq.  (A-27) 
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APPENDIX  II 

KIRCHHOFF  DIFFRACTION  BY  A  HALF-PLANE 


* 


I 


9 


A  plane  wave  field  u1  is 
x>0  as  indicated  in  Fig.  A-l . 


incident  on  the  half-plane  at  y=Q  for 
In  this  appendix,  the  expression  for  the 


Fig.  A-l.  Half-plane  geometry.  i 

field  diffracted  by  the  edge  (at  0)  of  the  half-pl ane  is  derived  by  using 
a  Kirchhoff  approximation.  Only  the  significant  details  of  the  analysis 
are  indicated.  It  may  be  verified  that  an  application  of  the  2-D  Green's 
theorem  to  the  region  surrounding  the  half-pl ane  leads  to  the  following 
integral  equation  for  the  electric  surface  current  density,  Jg  induced  on 
the  half  plane  when  a  plane  wave  consisting  of  a  z-directed  electric  field 
(TMZ  case)  is  incident  on  it[17]: 

m 

(A-28)  u(P)  =  -  !  G  (x,yjx\0)  J  (x',0)  dx' 
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where  GQ  is  the  2-D  free  space  Green’s  function,  and 


i*wJs(x',0)  =  |^t(x  '0+)  -  |jr(x  '°")  , 


(Note:  surface  current  density  J$  is  z-directed.)  u(x,y)  represents  the 
scattered  z-directed  electric  field  surrounding  the  half-plane.  The  geo¬ 
metrical  optics  approximation  allows  one  to  let 
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Employing  this  approximation  in  Eq.  (A-28)  gives: 
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(Note:  k  =  =  free  space  wave  number,  and  1Q  =  free  space  wave 

impedance  =  /p/e. )  The  representation  for  GQ  chosen  here  is  different 
from  that  used  in  Reference  [17]  for  solving  the  same  problem  (however, 
the  various  alternative  representations  for  Gq  may  be  obtained  via  the 
use  of  the  techniques  for  construction  of  Green's  functions  which  are  also 
given  in  Reference  [17]);  hence,  the  method  for  solving  Eq.  (A-29)  differs 
from  that  in  Reference  [17].  In  the  present  analysis,  the  following 
expression  is  employed  for  GQ: 
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Incorporating  Eq.  (A-30)  in  Eq.  (A-29)  gives: 
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in  Eq.  (A-31)  gives  the  result 
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(Here  k  =  k^  -  ik2  where  the  loss  k?  <<  k^  and  k-|  >  0.  The  loss  is  added 

2 

for  ease  of  analysis;  k2  may  be  set  equal  to  zero  in  the  final  solution.) 
Thus,  Eq.  (A-31)  becomes: 
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The  integral  in  Eq.  (A-32)  may  be  evaluated  by  the  steepest  descent  method 

which  includes  the  Pauli -Clenmow  modification[10]  (for  a  pole  near  a  saddle 

point)  after  employing  the  polar  transformations  k  =  kcosa,  and  k  =  +ksina. 

x  y 

The  integral  in  Eq.  (A-32)  is  thus  transformed  into  a  contour  integral  in 
the  complex  a-plane  as: 
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The  contour  of  integration  C  for  Eq.  (A-33)  is  shown  in  Fig.  A-2.  Since 
=  ir/2  in  the  present  problem,  Eq.  (A-33)  becomes  for  0<$<ir  (the  results 

A  A 

for  tt<4><2-it  are  symmetric  about  <|>=it): 
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The  integration  contour  C  may  be  deformed  onto  the  steepest  descent 
path  (SDP)  through  the  saddle  point  at  a=4»=ots  as  indicated  in  Fig.  A-2. 


Fig.  A-2.  Contour  of  integration  in  the  complex  a-plane. 


The  reflected  wave  pole  is  crossed  if  as<it/2  and  its  contribution  (residue) 
must  be  included.  The  saddle  point  contribution  to  SDP  gives  the  desired 
diffracted  field;  without  giving  details  of  the  SDP  evaluation  (which  are 
available  for  the  wedge  diffraction  problem  in  Reference  [10]),  the  result 
for  the  diffracted  field  u^  is: 
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where  F (X )  is  defined  in  Eq.  (23)  of  Section  II;  here  X  is  given  by: 
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Converting  the  result  of  Eq.  (A-35)  (which  is  valid  for  the  coordinate 
system  of  Fig.  A-l )  to  a  result  that  is  valid  in  the  coordinate  system  of 
Fig.  A-3  in  which  4>=4>-  and  $'=<)>'-  ^  =  0  (for  =  j) ,  one  obtains: 


Fig.  A-3.  Transformed  coordinate  system  for  the  half-plane  problem. 
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in  which,  the  Kirchhoff  diffraction  coefficient,  Dk(4>)  is 


(A38a) 

where 
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The  result  in  Eq.  (A-38a)  checks  with  that  given  in  Eq.  (22)  of  Section  II. 
The  ud(P)  for  the  TEZ  case  is  assumed  to  be  the  same  as  in  (A-37);  this 
assumption  is  reasonable  within  the  transition  region  near  4>=tt,  where  the 
Kirchhoff  result  is  relatively  insensitive  to  polarization  effects. 
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